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Abstract 

We report on the accuracy of the measurement of the two photon decay width of 
a Higgs boson with the mass of 120 GeV at the TESLA Photon Collider, assuming 
a 77 integrated luminosity of 80 fb _1 in the hard part of the spectrum. The QCD 
radiative corrections for the quark pair background processes are taken into account 
by a reweighting procedure. We found that the product r(H — ► 77) x BR(H — > bb) 
can be measured with the statistical error of 1.8% in one year of run. 



Introduction 



The search for the Higgs boson to explore the mechanism of symmetry breaking is one of the 
most important tasks of the experimental program at the future colliders. While it can only 
be produced in association with another particle at an e + e~ collider, the Higgs boson can be 
produced singly in the s-channel of the colliding photons at a photon collider. If the mass of 
the Higgs boson is light, as indicated by the current electroweak measurements, then it will 
be found by the time a 77 collider is constructed. The aim of such collider will be then a 
precise measurement of Higgs properties. Among these the measurement of the two photon 
decay width of the Higgs boson is particularly important. The Higgs boson can be produced in 
77 collisions via a one-loop diagram, in which all charged particles whose masses derive from 
the Higgs mechanism can appear inside the loop. In the Standard Model (SM), the dominant 
contributions come from the top and W boson loops. A deviation of the two photon width from 
the SM prediction indicates additional contribution from unknown particles, and is a signature 
of physics beyond the SM. For example, the minimal extension of the SM predicts the ratio of 
the two photon width T(h — ► 77, MSSM)/ T(H — > 77, SM) < 1.2 £Q for a Higgs boson with a 
mass of 120 GeV, assuming a supersymmetry scale of 1 TeV and the chargino mass parameters 
M and fi of 300 and 100 GeV, respectively. 

At a 77 collider, the two photon decay width is obtained by measuring the Higgs formation 
cross section in the reaction 77 — ► H — ► X, where X is the detected final state. The number of 
detected events is proportional to the product T(H — > 77)xBR(H — > X). Measuring the forma- 
tion cross section will determine this product. An independent measurement of the branching 
ratio BR(H — > X) at e + e~ collider allows a determination of the 77 partial width. 

In this paper we study the accuracy of the measurement of the two photon decay width for 
a Higgs boson with the mass of 120 GeV. The study is performed for the photon collider option 
of the TESLA Linear Collider operated at centre-of-mass energy of 210 GeV. The assumed 
data volume corresponds to a total 77 integrated luminosity of 400 fb _1 which amounts to 
about £ 77 (-y/i > 80GeV) = 80 fb _1 in the hard part of the spectrum. 

The feasibility of the measurement of the two photon decay width of the Higgs boson in 
this mass region has also been reported by [3]. 



Simulation of the signal and background processes 

The cross section for the Higgs boson formation is given by a Breit-Wigner approximation 
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where M H is the Higgs boson mass, r(H — > 77) and r tot are the two photon and total decay 
width of the Higgs boson, Ai and A 2 are the initial photon helicities and a/s 77 is the 77 centre- 
of-mass energy. The initial photons should have equal helicities, so that J z = 0, in order to make 
a spin-0 resonance as it is the case of the Higgs boson. If polarised photon beams are used, the 
signal cross section is increased up to a factor of 2. The experimentally observed cross section 
is obtained by folding this basic cross section with the 77 collider luminosity distribution. 

A Higgs boson with the mass of 120 GeV is detected in the 77 —>■ H — > bb final state. The 
event rate is given by the formula: 
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where the conversion factor (he) 2 is 3.8937966- 10 11 fb GeV 2 . 

It depends strongly on the differential luminosity, 4^-1 m„. The shape of the luminosity 

a V s 77 

distribution depends on the electron and laser beam parameters. The electron and laser beam 
energy considered for this study are 105 GeV and 1 eV, respectively, resulting in the maximum 
photon energy of about 75 GeV, suitable to study a Higgs boson with the mass of 120 GeV. 
The combination of the polarisation of the laser and the electron are chosen such to make 
the generated photon spectrum peak at its maximum energy. The helicity combination of 
the two high energy photons is arranged such that J z = state is dominant. The luminosity 
distributions for the J z =0 and J z = 2 are provided by the CIRCE2 [3] program and shown 
in Figure Beamstrahlung, secondary collisions between electrons and laser beam photons, 
and other non-linear effects are taken into account when calculating the 77 luminosity. The 
resulting value of j^" 7 |m h is 1-6 fb _1 /GeV. 

The branching ratios BR(H — > 77), BR(H — > bb) and the total width are taken to be 0.22%, 
68% and 4 MeV, respectively, and are calculated with the HDECAY [3] program including QCD 
radiative corrections. A signal rate of about 20000 events per year can be achieved under these 
conditions. 

The main background processes are the direct continuum 77 — > bb and 77 — > cc production. 
Due to helicity conservation, the continuum background production proceeds mainly through 
the states of opposite photon helicities, making the states J z = 2. Choosing equal helicity 
photon polarisations the cross section of the continuum background is suppressed by a factor 
Mq/s 77 , with M q being the quark mass. Unfortunately, this suppression does not apply to 
the process 77 — > qqg, because after the gluon radiation the qq system is not necessarily in a 
J z = state. The surviving background is large and overwhelms the signal. 

The signal 77 — > H — > bb process, as well as the background 77 — > qq(g) events are gener- 
ated with the PYTHIA jH] program. Parton evolution and hadronization are simulated using 
the parton shower and the string fragmentation models. A convolution with the luminosity 
distribution is performed, and kinematic cuts of | cos6>| < 0.9 and ->/s greater than 80 GeV 
are imposed during the event generation, 9 denoting the polar angle of the produced quarks. 
The QCD corrections to 77 — » qq have to be taken into account for a reliable background 
estimation. For this purpose we use the total cross sections shown in Figure El which include 
the soft and hard gluon emission [Jj, virtual corrections [Zj and non-Sudakov form factors [B] 
to weight the event. The event topology is given by the parton shower model. 

The effective cross sections, the number of expected and generated signal and background 
events are presented in Table [T] 

To simulate the response of the apparatus we use SIMDET , a fast Monte Carlo program 
which smears the kinematics of the particles in the final state according to the detector reso- 
lutions. The detector used in the simulation follows the proposal presented in the Technical 
Design Report of the TESLA e + e~ Linear Collider ^Hj ■ However, the changes in the detector 
design for 77 collisions will affect only the forward region ^1] which is not important for this 
analysis. 

Event selection 

The analysis aims to select events with two or three multiplicity jets from the Higgs boson 
decay. Two of these jets contain bottom quarks. The invariant mass of the jets has to be 
consistent with the Higgs mass. 
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Jets are reconstructed using the DURHAM clustering scheme with the parameter 

y cu t=0.02. 

We select events with a visible energy larger than 95 GeV and the energy imbalance longitu- 
dinal to the beam direction below 10% of the visible energy. FiguresE^b shows the distribution 
of the cosine of the thrust angle for the signal and background events. The s-channel signal 
process has an isotropic angular distribution, while the t-channel background processes are 
forward peaked. We require the absolute value of this quantity to be below 0.7. 

The cross section for the continuum production of the charm quark is 16 times larger 
than for bottom quarks, therefore b-quark tagging is crucial in this analysis. The b-tagging 
algorithm relies on a neural network with 12 inputs and 3 output nodes, described in Ref. |T3J. 
The most important inputs comprise the impact parameter joint probability tag introduced 
by ALEPH |14j . the pt corrected vertex invariant mass obtained with the ZVTOP algorithm 
written for the SLD experiment ^H] and a one-prong charm tag using the largest and second 
largest track impact parameter significances in r - <ft and r - z, in jets where ZVTOP found 
only the primary vertex. The spectrum of the neural network output NN out to discriminate 
b-quark jets from u-, d-, s- and c-quark jets is presented in Figure 13^,. The performance of the 
neural network b-tag in Z° decays is shown in Figure Eb- We require at least 2 vertices to be 
found and ask that the NN ont to be greater than 0.95. The b-tagging efficiency is 70% and 
corresponds to a purity of 98%. 

The efficiency of this selection is found to be 36%. 

Discussion of the QCD 77 — ► hadrons background 

A background not discussed yet is the process 77 — s>hadrons, where photons interact hadron- 
ically at the interaction point. With a total hadronic cross section in 77 collisions of 400 nb, 
we estimate about 1.4 such events per bunch crossing. These events obscure the interesting 
physics processes described in the previous sections. For this reason this class of events needs 
to be included in the PYTHIA simulation for overlap. The HADES JHI] program is used for 
this purpose. A large fraction of this background is distributed at small angles and we believe 
that it can be reduced cutting on the polar angle of the tracks [T7]. A careful study of this 
hadronic background is currently being performed. 



Results 

The reconstructed invariant mass for the selected signal and background events is shown in 
Figure 03 To enhance the signal a cut on the invariant mass is tuned such that the statistical 
significance of the signal over background is maximised. Events in the mass region of 106 GeV 
< Mjj < 126 GeV are selected. The number of estimated signal and background events in this 
window are 6018 and 7111, respectively. 

The two photon decay width of the Higgs boson is proportional to the event rates of the Higgs 
signal. The statistical error of the number of signal events, VN bs/(N obs — iV bkg ), corresponds 
to the statistical error of this measurement. Here iV b s is the number of observed events, while 
iVbkg is the number of expected background events. 

We obtain 

A[r(H^77)xBR(H^bb)] r- 
[T(H - 77) x BR(H - bb)] " VA W(^ob s - ~ 1-9%- 



4 



We can improve this number by 0.1% doing a fit with the function 

, _ y iifjNH + kl-^l 

X — ^ j ) 

di 

with fi — ^ being the fraction of signal events, hi and di the background and observed number 
of events in each bin of the invariant mass distributins and Nr a fit parameter. 

We conclude that for a Higgs boson with a mass Mh=120 GeV we can measure the prod- 
uct T(H — > 77) x BR(H — > bb) with an accuracy of 1.8% using an integrated luminosity 
corresponding to one year of data taking at the TESLA Photon Collider. 
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Table 1: Cross sections and the number of expected and generated events for the signal and 
background processes. 
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Figure 1: Luminosity spectra for the spin of the two colliding photons J z = and J z = 2. 
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Figure 2: Cross sections for the 77 
corrections for J-, = and J z = 2. 



bb(g) and 77 — > cc(g) processes including QCD radiative 
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Figure 3: Distributions for the cosine of the thrust angle for signal a) and background b) events. 
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Figure 4: a) Distribution of the neural network output and b) efficiency versus purity curve for 
the neural network based b-tag in Z° decays. 
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Figure 5: Distribution of the reconstructed invariant mass for the signal and background events. 
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